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Abstract 

For a dense stellar matter, which is electrically neutral and in beta equilibrium, the 
electron chemical potential, u e , will depend nontrivially on baryonic matter density. It is 
generally expected that as density increases, the electron chemical potential will increase and 
new degrees of freedom will emerge as \x e becomes comparable to their energy scales. As- 
suming the electrical neutrality and beta equilibrium for the stellar matter, we have studied 
how the density dependence of lepton chemical potentials varies for different models of nu- 
clear interactions that are constrained by experiments up to nuclear matter density, no, but 
extrapolate differently(unconstrained) beyond uq and calculated the relative abundances of 
nucleons (neutron and proton) and leptons (electron and muon) and their density dependen- 
cies. We find that the density dependence of the electron chemical potential is strongly 
dependent on the structure of the nuclear symmetry energy relevant to softness/stfness of 
the nuclear matter EOS that measures the energy relevant to the neutron-proton asymme- 
try. As a consequence, the relative abundances of neutrons, protons, electrons, and muons 
as well as the kaon condensation are strongly dependent on the nuclear symmetry energy. 
An intriguing result in our finding is that contrary to the accepted lore, kaon condensation 
in neutron star matter, which is considered to be the first phase transitions beyond no and 
plays a crucial role in certain scenarios of compact-star formation, is not directly tied to the 
softness or stiffness of the EOS beyond no- This point is illustrated with a "super-soft" EOS 
that is fit to the ^ ratio data of GSI which excludes kaon condensation at any density. 



1 Introduction 



Recent interest on the astrophysical compact objects, i.e., neutron stars and black holes, and 
their evolutions (cooling and collapse for example) brings out a very challenging question 
as to what is the relevant equation of state(EOS) of the stellar matter, deep inside which 
the density is expected to be much higher than the normal nuclear density, no- There have 
been lots of development in constructing the equations of state for nuclear matter, which are 
model dependent but supposed to be consistent with the experiments up to no [I]. However, 
the simple extrapolations of equations of states beyond uq lead to different and strongly 
model-dependent results for the macroscopic variables: chemical potentials of constituents, 
energy density, pressure, etc.. And it becomes a very important subject to determine or 
predict EOS far beyond normal nuclear density. 

The difference of stellar matter from the nuclear matter is that there are additional 
constraints for a stellar matter, which are the conditions of electrical neutrality and beta 
equilibrium. These constraints naturally lead to the emergence of leptons (electron and 
muon) as well as protons and the lepton chemical potential is expected to be increasing 
with the density. It opens new degrees of freedom for the stellar matter when the scale for 
the new degrees of freedom is comparable to the lepton chemical potential. 

For an ideal noninteracting fermi gas of neutron and proton the lepton chemical 
potential is determined by the Fermi-Dirac statistics only and is known to be increasing 
slowly with the nucleon density. But the nucleon system is very strongly interacting and 
the isospin asymmetric nuclear interaction, which is known as symmetry energy of nuclear 
matter, cannot be simply ignored in the stellar matter. The symmetry energy is defined 
by the sum of the kinetic contribution and asymmetric nuclear interactions, which are 
required to have neutron-proton asymmetry. The density dependencies of the symmetry 
energy have been discussed in the various theoretical models [lj [2 J |3j and also have been used 
for the phenomenological applications including heavy ion experiment[4j and astrophysics. 
Among the interesting issues, where the importance of symmetry energy in astrophysical 
phenomena has been well addressed, are the size of neutron stars[5][6][7][8], the onset of 
kaon condensation [9j [TU\, the cooling of proto-neutron star[TT], gravitational waves [T2] and 
the instability of neutron stars |13j. 

In this work we focus on the density dependent change of relative abundances of 
the stellar matter with particular emphasis on role of the nuclear symmetry energv[T3]|15|. 
which we try to elaborate more transparently using the various models of symmetry energy. 
As noted above, as density is increasing, we are expecting the emergence of new degrees 
of freedom very essential in formulating the corresponding EOS, which is only possible 
when we know the abundances of the constituents. Also we could get a more transparent 
understanding of EOS in terms of its constituents. 

The phenomenological models of the symmetry energy show similar density depen- 
dencies up to near the normal nuclear density but their density dependencies diverge quite 
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widely beyond the nuclear density from model to model(for a recent review see Ref. [lj). 
Therefore any simple extrapolation far beyond normal density should result in quite dif- 
ferent conclusions, which might be the case for investigating stellar matter with higher 
density. In this work we calculate the density dependence of the abundances using a few 
selective models of nuclear symmetry energy, which have been used extensively to calculate 
neutron star EOS in detail[lj|7j|9j[10j|16j. We demonstrate how the lepton chemical po- 
tentials depends on the nuclear density and also how the results depend on models of the 
nuclear symmetry energy. It is worth noting that ir~ /tt + ratio in heavy ion collisions has 
been known to be sensitive on the symmetry energy. Recent studies by Xiao et al. [T7] of 
FOPI data at SIS/GSI seem to show that the super-soft EOS that gives the x=l curve in 
|17j is the only EOS so far available that fits the pion data of GSI. However the symmetry 
energy factor with x=l drops to zero before reaching higher densities than ~ 3no, which 
on the other hand cover the relevant range of density in a stellar matter. Hence the beta 
equilibrium in stellar matter might not be expected nor the kaon condensation as discussed 
below. 

When the density is increasing in stellar matter new degrees of freedom, excited 
baryons and mesons including hyperons, are supposed to emerge. Kaon is among the 
possible new degrees of freedoms after light leptons, when the medium dependent energy of 
kaon becomes comparable to the lepton chemical potential. Employing the simple formulae 
for the inverse propagator of kaon in medium [18], the threshold densities for the kaon 
condensation are calculated using the model-dependent symmetry energies. It is a very 
interesting issue since Bethe-Brown argument [19J for the maximum mass of neutron star , 
Mfflg ~ 1.5Mq, depends on the onset of the kaon condensation at the density of ~ 3reo- 

In this work, the density dependence of the relative abundances of stellar constituents ( 
neutron, proton, electron and muon)are calculated using several models of density depen- 
dent nuclear symmetry energy, since it is directly related to the relative abundances of 
constituents. For the possible implication of the kaon condensation, the threshold densities 
of kaon condensation are calculated by comparing the lepton chemical potentials with the 
zero of the kaon inverse propagator. 

Let us start with a simple matter consisting of pure neutron gas. The temperature 
is taken to be zero throughout this work. The energy density and pressure of a simple free 
neutron gas[20)(We put H = c = 1) are given by : 
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where pp is the Fermi momentum defined by the neutron number density 
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The Fermi momentum in unit of normal nuclear density , no = 0.16fm 3 , and the kinetic 
energy(Fermi level) is given by 



/ \ 1/3 2 / \ 2/3 

PF = 336MeV I — , E k F in = = 60MeV — . (1.4) 

\n J ' 2m N \noJ 

The chemical potential of neutron, [i n , is the same as the Fermi energy Ep : 

fi n = E F = (p 2 F + mlf 12 = m (1 + x 2 ) 1/2 ■, (1-5) 
where x is a dimensionless parameter defined by 

x = — . (1.6) 

m 

It is interesting to note that there is isospin symmetry in nature. A proton is the 
isospin partner of the neutron, which consists of an iso-doublet with almost the same mass, 
~ 940MeV. Hence if there is any channel for the neutron to be converted to proton, the 
energy of the system of neutron only can be lowered. The conversion depends on the 
dynamics, which is supposed to be the /3-decay, 

n—>p + e + b> e (1.7) 

to reach /3-equilibrium eventually for which 

Hn = Hp + He + Hv (1-8) 

Hereafter we will assume all neutrinos are emitted out of star, then we take the neutrino 
chemical potential as fi u = 0, fi n = u p + fi e . 

If the energy only is concerned, then the minimum can be achieved when u n = u p 
without electron, \i e = . However, a stellar matter is believed to wind up as a charge 
neutral object. Practically during the period of reaching the (3 equilibrium, the stellar 
matter reacts to satisfy the neutrality condition, 

n-p = n e or m p x p = m e x e , (1-9) 

where x p and x e are the dimensionless Fermi momentum of proton and of electron, respec- 
tively. This implies the number of protons should be balanced by the number of electrons, 
which means the energy minimum can be reached with the constraint of neutrality for a 
system of free neutron, proton and electron(NPE gas) in a (3 equilibrium. It is one of 
the essential differences from nuclear matter, where the charge neutrality is not a physical 
constraint, for example, for heavy ions. 

As the density increases such that the electron chemical potential is comparable 
to muon mass, m M , it is energetically favorable for the electrons to convert to muons as 
e — > fx + v e + Dp. Then the chemical equilibrium can be accomplished as 

- «p = Me = AV, ( L1 °) 
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Figure 1: (a)The relative abundances of particles and (b)/x e vs. u. Solid line and dashed 
line respectively refer to NPE/x gas and NPE gas. 



with the neutrality condition, 

n p = n e + n ll . (1-H) 

At the muon threshold, the muon density is zero, = or = and = = 106MeV 
and the Fermi momentum of electron at threshold is determined by the masses of muon 
and electron as 

1/2 

1 



X 



e thres 



, 2 

nifj, 



m, 



207. (1.12) 



The nucleon density at threshold can be determined by Eq. (jl.lip . n p = n, 



71, 



nL-thres = 2.91n , — L-thres = 0-011. (1.13) 
n 

Beyond the muon threshold the constituents of the stellar matter becomes neutron, proton, 
electron and muon(NPE^ gas). 

The relative abundances of neutron, proton, electron and muons are shown in Fig. 
[TJa). The relative neutron abundance is almost not changing but the relative abundances 
of proton are increasing substantially beyond the muon threshold as expected. The proton 
fraction, n p /n, increases with density as well as for the electron, n e /n. 
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The electron chemical potential increases with density as shown in Fig. Q3h). In 
the NPE gas, the number of electrons is equal to that of protons by the charge neutrality 
condition, Eq. (jl.9p . Until the muon threshold density(~ 3no), the increment of the electron 
chemical potential follows the property of the NPE gas. Then, beyond the threshold density, 
the electron chemical potential grows with density by the property of the NPE// gas. In 
the NPE/i gas electron shares its number with muon such that the number of electron is 
reduced by the presence of muon. 



2 Relative abundances with nuclear symmetry energy 

So far, in calculating the chemical potentials, the nucleons are considered as free particles, 

- ^ee = mn(1 + ^2)1/2 _ ^ + ^2)1/^ ^ 

which can not be realistic due to the nuclear interactions. 

It is useful to write the energy per particle in the nuclear matter [2T] as 



o / \ 2/3 

E(n,N p ) ~m N + -E° F f — + S(n)(l - 2N p f + V(n), (2.2) 
5 V n o/ 



where the last term does not depend separately on neutron or proton number density but 
only on total number density. In the second term, E F = — — ^ s ^ ne F errrn energy at, 
n = no/2. The third term is the symmetry energy for nuclear matter 

E sym = S(n)(l-2N p ) 2 , (2.3) 

and the nuclear symmetry energy density, 

e sy m = nS(n) (1 - 2iV p ) 2 , (2.4) 



where the symmetry energy factor, S(n), is a model dependent function of density [16J [21 J [22J 
[23] . which is manufactured to be consistent with nuclear matter data up to the nuclear 
matter density, no- 

Now we can calculate the contributions of the symmetry energy density, e sym , to the 
chemical potentials of proton and neutron to get 

Mn - u p = fir - n;*™ = 4(1- 2N p ) S(n). (2.5) 

We can see that the larger S(n), the lager the chemical potential difference. 

For an ideal free gas, the chemical potential difference, in the non-relativistic limit of 
Eq. ([13]) is given bv[24] 



- 4 ree = (3vr2n)2/3 i (1 " Np)2/3 ~ n2 A ' (2 - 6) 
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which can be approximated in the well known form|17|. 

M /ree _ ^ree = 4 (1 _ 2A y Sfree (n), (2.7) 

where 

S fr Un) = (2^-l)\E%(^'\ (2.8) 

There are a number of parameterizations for the symmetry energy factor, S(n), which take 
account of the free kinetic contribution as well as the potential energy contribution. One 
of the parametrization used in Refs. [16] [21] is 



S F {n) = (2^-l)\E% 
5 



\ 2 / 3 

- -Fin) 
no J 



+ S F(n), (2. 



where F(n) is a parameter for the potential contribution to the symmetry energy, which 
satisfies F(0) = and F(tiq) = 1. Hereafter we set F(n) simply as F(n) = ti/uq. In the 
last term of Eq. (|2.9|) . So denotes the bulk symmetry energy parameter, So ^ 30MeV. 
One of the different forms of S(n) suggested in [22] is 

o / \ 2/3 / \ B(a) 

S a = ^-l)p° F (^) + A(a)- + [18.6 -A(a)} - , (2.10) 
5 \n J n \n J 

where a is a free parameter which determines the bulk property of nuclear matter. In this 
work, we take a = 1 , which reproduces ir + /n~ ratio in heavy ion collision |17j. 

A(a = 1) ~ 107MeV and B(a = 1) ~ 1.25 (2.11) 

(for other parameter sets, see Ref. |22|). In Eq. (|2.10|) . the first term denotes the kinetic 
contribution and last terms refer to the potential energy contribution of nuclear matter as 
in Sp. 

There is also a different type of symmetry energy factor, in which different parametriza- 
tion scheme is adopted as [23 1, 

S 3 (n)^S* +Lp + ±Kp 2 , (2.12) 

where 

n — n* . 

>S?-> (2 ' 13) 

with n* = 0.148fm~ 3 and £3 represents either SVsu or 5nl3 in Ref. [23]. The bulk pa- 
rameters, Sq, L and K in Eq. (|2.12p consistent with the nuclear matter are listed in Table 
1. 

For comparison, we plot Sf, S a= i, Sfsu and 5nl3 in Fig. [2]. As one can see these 
symmetry energy factors are not much different from each other up to uq. But beyond no, 
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Model 


So 


L 


K 


FSU 


32.59 


60.5 


-51.3 


NL3 


37.29 


118.2 


100.9 



Table 1: Bulk parameters for the symmetry energy factor in Ref. |23| . 
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Figure 2: The density dependencies of 5(n)'s for different models. 

the difference becomes significant from model to model as shown in Fig. [2j In general the 
symmetry energy factors are increasing with the density up to n ~ 5riQ. However S a= i has 
peak near n ~ no and drops afterwards to zero at n ~ 3no- 

The electron chemical potential and its density dependence in /3-equilibrium can be 
calculated using Eqs. (jl.lOp and (|2.5p for different symmetry energy factors, as shown in 
Fig. [3j The electron chemical potential for the free neutron gas is lower than with other 
symmetry energy factors, except S a , beyond n ~ 2uq. This is the main reason why it is 
difficult to excite new degrees of freedom in the lower density for the case of the free nucleon 
gas. For S a =i, the density dependence of the electron chemical potential is quite different 
such that it drops to zero n ~ 3no- It implies that /3-equilibrium may not be reached for 
n > 3no with the symmetry energy factor of S a= \. 

When the electron chemical potential increases such that it becomes comparable to 
the rest mass of muon, muons become constituents of stellar matter in addition to neutrons, 
protons and electrons. At the muon threshold with n M = the electron chemical potential 
is just the mass of muon as discussed in previous section. Then with given proton and 
electron densities, we can determine the nucleon density, «L-thres> a t muon threshold from 
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Figure 3: The electron chemical potentials, \x e due to different models for sNPE gas. 





Free 


S F 


5 FSU 


%L3 5 '«=i 


n ^-thres 


2.91 


0.94 


0.76 


0.73 0.73 



Table 2: The muon threshold densities obtained using different models 
for S(n), in unit of no- 

the chemical equilibrium condition, given by 

A{l-2N p )S{n) = m e (l + xl) 1 ' 2 . (2.14) 

In the previous section, we get n M _ t hres = 2.91no in Eq. (|1.13p for free nuclear matter. 
With the symmetry energy factor, Sp{n), Eq. (|2.14j) gives 

^U-thres = 0.94n , (2.15) 

which is much lower than for the free case. Threshold densities for other symmetry energy 
factors can be easily guessed in Fig. [3] and the corresponding threshold densities are calcu- 
lated using Eq. (|2.14p in Table 2. One can note that the muon threshold density becomes 
about three to four times lower than for the free nuclear matter. 

At higher density beyond the muon threshold, the stellar matter consists of neutron, 
proton, electron and muon(sNPE/x gas). And the chemical equilibrium condition and charge 
neutrality condition become 

4(1 - 2N p )S(n) = m e (l + x2) 1 /2 = m/t ( 1 + x 2 ) i/ 2) (216) 

n p = n e + n^. ( 2 -17) 
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Figure 4: The electron chemical potentials for sNPE/z gases with different nuclear symmetry 
energy models. 



Effectively these are three equations to be solved for four unknowns, x e , i„, x n and x p . Then 
for a given electron chemical potential, x e , one can solve these equations completely to de- 
termine the chemical potentials and the relative abundances of the constituents, n e ,n^,n n , 
and n p . The electron chemical potentials(same as the muon chemical potentials in beta 
equilibrium) beyond the muon threshold densities(see Table 2) are calculated for different 
models of the nuclear symmetry energy as plotted in Fig. SJ The density dependence of 
electron chemical potential for Sfsu shows a quite different behavior from others. One 
can observe in Fig. [5] that the proton and muon fractions are increasing with the density 
much faster than in free case. For Sfsu, the fractions of proton and muon are increasing 
with the density up to ~ 4no but decreasing and even becomes lower than free case for 
the higher density. The S a case shows a similar behavior with Sfsu, but the lepton and 
proton fractions are decreasing more quickly than Sfsu an d the muon and proton fractions, 
beyond the muon threshold density, are much lower than in other models. 

When the leptonic chemical potentials increase high enough, it is natural to ask what 
kind of new degrees of freedom can be driven to evolve in a stellar matter. They could 
be excited baryons including hyperons as well as strange mesons. One of the interesting 
possibilities is the s-wave kaon condensation [25] driven by the electron chemical potential. 
The simplest way of demonstrating the possibility of kaon condensation is by comparing 
the electron chemical potential with the kaon chemical potential at the threshold, which 
determines the kaon threshold density. At the threshold the kaon chemical potential is given 
by the zero of the inverse propagator, D^_, with zero momentum[18] which is given by 
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Figure 5: The relative abundances of constituent particles for sNPE/x gases with different 
nuclear symmetry energies. 



D K- = U K ~ m K + -pi 71 ™/ 2 + n p)uK + — p~ n , (2.18) 

where / is the pion decay constant, / = 93MeV. The possible range of £^at[18J[26J[27J[28J 
is estimated to be 

T, KN = 200 - 400MeV. (2.19) 

As one of the example, we take T>kn = 400MeV in this work[9j. It is found that up to 
n ~ 5no the dependence of ojk (n) on the symmetry energy models is not so significant and 
the kaon threshold densities are estimated around n ~ 3no as shown in Table 3 and in Fig. 
EJ The threshold density is lower ~ 2.5no with Snl3 but higher for free case, ~ 3.5no- One 
should note that lor obtained using Eq. (12. 18ft is found to be always larger than electron 
chemical potential such that one cannot expect a kaon threshold if S a= i is valid up to 
n ~ 3uq. 
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Figure 6: The electron chemical potentials and kaon chemical potentials for different nuclear 
symmetry energy models. 



3 Discussion 

In this work we consider the effect of the nuclear symmetry energy on the relative abun- 
dances of particles, neutron, proton, electron and muon, with the nucleon density inside a 
stellar matter supposed to be electrically neutral and in beta equilibrium. We determine 
the electron and muon threshold densities. The muon threshold density reduces substan- 
tially with symmetry energy from that of free nucleon gas. It is observed that the relative 
abundances with symmetry energy diverge significantly from model to model for higher 
density n > no- We also estimate the kaon condensation threshold densities using a number 
of different models for the symmetry energy. The kaon threshold estimated with symmetry 





Free 


Sf 


Sfsu 


SnL3 


n [n ] 


3.53 


2.76 


2.89 


2.45 


H [MeV] 


120.4 


206.3 


192.5 


234.5 



Table 3: The kaon threshold densities estimated using different models 
for S(n) and the electron chemical potentials at the threshold densities. 
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energy considered is found to be lower than the free case. However, for the super-soft EOS, 
in which the symmetry energy factor S a= i drops to zero before reaching higher density 
than ~ 3uq, beta equilibrium in stellar matter cannot be expected. Moreover u>k obtained 
using Eq. (|2.18p is found to be always larger than electron chemical potential such that one 
can not find kaon threshold for S a= \. This is an intriguing consequence of our finding that 
contrary to the accepted lore, kaon condensation in neutron-star matter which is considered 
to be the first phase transition as density increase beyond no and hence plays a crucial role 
in certain scenarios of compact-star formation is indifferent to the softness/stiffness of EOS 
beyond no- 

We can notice that the presently known models for the symmetry energy factor pro- 
vide quite different predictions for higher densities when the parametrization is straightfor- 
wardly extrapolated to higher densities beyond ~ 3no, which are however very important 
deep inside the stellar system. Hence the predictions on the stellar matter structure after 
integrating the Tolman-Oppenheimer-Volkov equation should depend strongly on the spe- 
cific form of symmetry energy[6j[7j[10j[29j. Hence it is very important to find the symmetry 
energy with correct density dependency for a stellar matter and we are looking forward to 
get more experimental information on symmetry energy from the forthcoming experiment 
like FAIR/GSI. 

One of the interesting ideas developed recently in hadron physics is adopting the 
hidden local symmetry |30|. which opens up an idea of hadronic freedom[31J. However it 
is not well investigated how this idea can be employed in understanding the symmetry 
energy at higher density and its relation to the softness/stifness of EOS. Whether it simply 
converges to the free case for higher density than ~ 3uq or it is quite different matter from 
the free fermi gas and how the symmetry energy is related to the KN interaction which 
controls the kaon condensation at high density nuclear matter are among the very intriguing 
questions to be answered. Since it is the region so far unexplored by the experiment, it is 
very interesting to construct a form of symmetry energy compatible with the new idea of 
dense hadronic matter to investigate the implications on the stellar structure using TOV 
equation, which will be discussed in a separate work. 
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